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LAMINATION  RESIDUAL  STRESSES  IN  MULTILAYERED  FIBER  COMPOSITES 

by  Christos  C.  Chamis 
Lewis  Research  Center 


SUMMARY 


Residual  stresses  arising  from  the  lamination  fabrication  process  of  composites 
are  investigated  using  linear  laminate  theory.  The  resulting  equation  for  predicting 
residual  stresses  is  presented.  The  variables  that  influence  the  residual  stresses  are 
identified.  These  variables  are  constituent  materials  and  their  thermal  and  elastic 
properties,  average  fiber  diameter  and  number  of  fibers  per  end,  fiber  and  void  vol¬ 
ume  ratios;  ply  stacking  sequence  and  orientation,  and  processing  temperature.  The 
effects  of  these  variables  on  the  residual  stresses  of  several  composite  systems  are 
examined.  The  fiber  composites  considered  are  boron/aluminum,  boron/epoxy, 

S- glass/epoxy,  Thornel-50  graphite/epoxy,  and  Modmor-I  graphite/polyimide. 

It  was  found  that  the  residual  transverse  tensile  and  in-plane  shear  stresses  can 
reach  magnitudes  comparable  to  the  corresponding  ply  strengths  and  can  cause  transply 
cracking.  The  ply  residual  stresses  are  very  sensitive  to  constituent  material  proper¬ 
ties,  ply  stacking  sequence  and  orientation,  fiber  content,  and  processing  temperature. 
They  are  relatively  insensitive  to  the  void  content.  The  residual  stresses  are  relatively 
small  in  composites  where  the  fiber  content  or  the  processing  temperature  vary  only 
moderately  from  ply  to  ply.  Ply  residual  stresses  increase  with  increasing  matrix 
modulus.  They  increase  linearly  with  increasing  matrix  thermal  coefficient  of  expan¬ 
sion.  It  is  possible  to  orient  plies  in  simulated  components  so  that  some  plies  have 
tensile  longitudinal  residual  stress.  Residual  stresses  can  cause  interply  delamination 
due  to  the  relative  rotation  of  adjacent  plies. 

Residual  stress  transply  cracking,  in  general,  can  be  prevented  by  increasing  the 
ply  transverse  tensile  and  shear  strengths,  decreasing  the  processing  temperature,  or 
decreasing  the  difference  between  the  fiber  and  matrix  thermal  coefficients  of  expansion. 
In  specific  designs  transply  cracking  can  be  prevented  by  suitable  ply  orientation,  in¬ 
crease  in  fiber  content,  use  of  transitional  plies,  and  reduction  of  the  number  of  off- 
axis  plies. 

Experimental  work  to  fully  assess  the  residual  stress  effect  on  the  structural  re¬ 
sponse  of  the  composite  is  suggested. 


INTRODUCTION 


Multilayered  fiber  composites  are  readily  adaptable  to  simultaneous  component  and 
material  designs.  This  approach  allov/s  the  designer  to  take  full  advantage  of  the  high 
stiffness  and  strength-to-density  ratios  of  fiber  composites.  The  result  is  an  efficient 
use  of  fiber  composite  material  in  structural  components.  The  composite  is  usually 
designed  by  specifying  the  fiber  direction  in  each  ply,  the  number  of  plies,  and  the  fiber 
volume  ratio  necessary  to  meet  the  given  loading  requirements. 

However,  having  different  fiber  directions  in  different  plies  can  cause  transply 
cracks  such  as  have  been  observed  in  graphite -fiber  composites.  These  cracks  are 
caused  by  the  large  thermal  residual  stresses  present  in  the  composite.  The  residual 
stresses  result  from  the  fabrication  process  used  to  mxake  fiber  composites.  The  re¬ 
sidual  stresses  arising  from  the  fabrication  process  are  examined  quantitatively  herein 
using  linear  laminate  theory. 

Multilayered  fiber  com.posites  consist  of  several  unidirectional  layers  (plies) 
oriented  at  specified  angies  relative  to  some  coordinate  axes  system.  The  geometry  of 
such  a  composite  is  depicted  in  figure  1.  The  plies  consist  of  parallel  fibers  embedded 
in  a  matrix.  The  geometry  of  a  ply  is  shown  in  figure  2.  Glass  or  graphite  fiber  com¬ 
posites  will  have  many  filaments  through  the  ply  thickness. 

In  general,  composites  achieve  their  structural  integrity  as  follows:  Nonmetallic 
matrix  composites  are  cured  at  elevated  temperature  under  pressure;  metallic  matrix 
composites  are  either  hot  pressed  or  plasma  sprayed.  These  methods  induce  residual 
stresses  in  the  composite  when  it  is  at  a  temperature  different  from,  its  cure  or  proc¬ 
essing  temperature.  The  difference  between  the  cure  or  processing  and  the  use  tem¬ 
perature  will  be  referred  to  here  as  temperature  difference.  The  residual  stresses  can 
be  either  microresidual  or  macroresidual  stresses. 

The  microresidual  stresses  are  influenced  by  the  local  geometry  of  the  constituents 
within  the  ply.  They  arise  from  the  difference  in  the  thermal  coefficients  of  expansion 
of  the  constituents  and  the  temperature  difference.  The  matrix  phase  change  shrinkage 
for  resins  also  contributes  to  the  microresidual  stress.  Adams,  Doner,  and  Thomas 
(ref.  1)  and  Haener,  Ashbaugh,  Chia,  and  Feng  (ref.  2),  for  example,  treat  this  type  of 
residual  stress  analytically.  Marloff  and  Daniel  (ref  3)  and  Koufopoulos  and  Theocaris 
(ref.  4)  examined  microresidual  stress  experimentally.  The  ply  shown  in  figure  2  will 
have  microresidual  stresses. 

Macroresidual  stresses  (ref.  5)  are  integrated  averages  through  the  ply  thickness. 
These  stresses  have  an  effect  on  the  ply  similar  to  an  equivalent  thermal  or  mechanical 
load  applied  to  the  composite.  They  arise  from  the  difference  between  the  ply  longitu¬ 
dinal  and  transverse  thermal  coefficients  of  expansion.  Macroresidual  stresses  are 
present  in  those  composites  that  have  plies  oriented  at  different  angles  and  whose  use 
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temperature  differs  from  the  cure  or  processing  temperature.  These  stresses  will  be 
present  in  the  composite  depicted  in  figure  1. 

The  significance  of  the  lamination  residual  stress  was  realized  when  transply 
cracks  were  observed  in  0*^/90*^  graphite  composites.  Doner  and  Novak  (ref.  5)  pre¬ 
sented  photomicrographs  showing  transply  cracks.  DeCrescente  and  Novak  (ref.  6)  and 
Winters  (ref.  7)  have  observed  similar  cracks  in  other  graphite/epoxy  composites. 

Some  analytical  treatment  of  the  lamination  residual  stress  has  also  been  reported. 
Doner  and  Novak  (ref.  5)  considered  the  residual  stresses  in  0°/90°  composites. 
Edighoffer,  Ravenhall,  and  Juneau  (ref.  8)  examined  the  residual  stresses  in  a  boron/ 
epoxy  composite  with  laminations  ±6/90  where  the  ±0  are  considered  as  an  equivalent 
0°  ply.  This  author  (ref.  9)  used  linear  laminate  theory  to  predict  residual  stresses  in 
graphite  composites  without  restriction  to  ply  orientation.  The  lamination  theory  used 
in  reference  9  was  developed  in  reference  10.  It  was  implemented  in  references  11 
to  13  and  in  a  computer  code  (ref.  14). 

The  computer  code  of  reference  14  provides  a  method  for  micromechanics,  lami¬ 
nation,  and  stress  analyses  of  multilayered  fiber  composites.  The  inputs  to  this  code 
are  constituent  material  properties,  fiber  and  void  contents,  correlation  coefficients, 
composite  geometry,  temperature  difference  (use  temperature  minus  cure  tempera¬ 
ture),  and  composite  stress  resultants  or  displacements. 

The  computer  code  of  reference  14  permits  the  variation  of  several  basic  param¬ 
eters  in  fiber  composites.  It  is  used  herein  to  investigate  the  residual  stresses  in  lami¬ 
nated  fiber  composites.  The  following  systems  are  examined:  boron/aluminum,  boron/ 
epoxy,  S-glass/epoxy,  Thornel- 50  graphite/epoxy,  and  Modmor-I  graphite/polyimide. 
The  significant  parameters  in  these  systems  are  ply  orientation,  fiber  and  void  contents, 
matrix  modulus  of  elasticity,  matrix  and  thermal  coefficients  of  expansion,  and  the  tem¬ 
perature  difference.  To  the  author's  knowledge,  the  effects  of  some  of  these  param¬ 
eters  on  the  residual  stresses  are  investigated  for  the  first  time.  Also  investigated  for 
the  first  time  are  (1)  possible  ply  delamination  due  to  residual  stress,  and  (2)  residual 
stresses  in  composites  which  are  unsymmetric  with  respect  to  bending. 

The  ply  orientations  examined  are  shown  in  table  I  together  with  the  notation  adopted 
herein  for  convenience.  Theoretical  thermal  and  elastic  properties  for  the  composite 
systems  considered  are  listed  in  table  II.  These  properties  are  for  a  fiber  volume  ratio 
of  0.  55.  The  corresponding  theoretical  strength  properties  are  given  in  table  III.  The 
results  in  tables  n  and  III  are  those  predicted  by  the  computer  code  (ref.  14). 

The  significant  parameters  mentioned  previously  are  varied  in  the  computations. 

The  results  are  plotted  to  illustrate  the  dependence  of  the  ply  residual  stresses  on  these 
parameters  as  is  described  subsequently.  The  plots  are  grouped  by  composite  system. 
The  composite  systems  boron/aluminum,  boron/epoxy,  S-glass/epoxy,  and  Thornel-50 
graphite /epoxy  are  used  in  some  basic  ply  orientations.  Results  for  various  ply  orien- 
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tations  in  composites  unsymmetric  with  respect  to  bending  are  also  presented  for  these 
systems.  The  composite  system  Modmor-I  graphite/polyimide  is  used  in  simulated 
structural  components.  Some  results  are  presented  for  this  system  to  demonstrate 
possible  ply  de lamination  due  to  residual  stresses. 

The  results  presented  herein  are  for  three  temperature  differences:  900®  F  (773  K) 
for  boron/aluminum,  300®  F  (440  K)  for  the  boron/epoxy,  S- glass/e  poxy,  and 
Thornel-50  graphite/epoxy  systems,  and  600®  F  (606  K)  for  the  Modmor-I  graphite/ 
polyimide  system.  Residual  stresses  for  any  other  temperature  differences  are  ob¬ 
tained  from  the  plots  by  direct  proportion.  This  is  possible  because  a  linear  theory  is 
used  and  the  residual  stress  is  linearly  dependent  on  the  temperature  difference. 

Several  recommendations  are  made  for  minimizing  or  eliminating  the  residual 
stress.  These  recommendations  are  deduced  from  the  results  presented  and  the  equa¬ 
tions  used  to  predict  the  residual  stresses.  These  equations  are  presented  in  their 
final  form.  The  derivations  are  given  in  reference  10. 


SYMBOLS 

Acx  array  of  composite  extensional  stiffnesses  referred  to  composite  structural 

axes 

Ccx  array  of  composite  coupling  stiffnesses  referred  to  composite  structural  axes 

Dcx  array  of  composite  bending  stiffnesses  referred  to  composite  structural  axes 

array  of  strain-stress  relations  (elastic  constants  of  the  i*^  ply) 

E^ll  ply  longitudinal  modulus 

E^22  ply  transverse  modulus 

0^12  ply  (in-plane)  shear  modulus 

fiber  and  void  volume  ratios,  respectively 

^cATx  vector  of  unbalanced  thermal  moments  referred  to  composite  structural  axes 

^cATx  vector  of  unbalanced  thermal  forces  referred  to  composite  structural  axes 

array  of  transformation  coefficients  for  the  i**^  ply 
S^llT  ply  longitudinal  tensile  strength 

S^llC  ply  longitudinal  compressive  strength 

S^22t  transverse  tensile  strength 

S^22C  P^y  transverse  compressive  strength 
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^Z12S 

x,y,z 

^li 

1,2,3 

“Zi 

«Z11 

"Z22 

%2T 

^cox 


^Zi 

^mpT 

e 


«cx 

‘'Z12 


•^Zi 

‘^/i22 
{  } 

[  ] 


ply  in-plane  shear  strength 
temperature  difference  for  the  i^^"  ply 
structural  axes 

th 

distance  from  reference  plane  to  centroid  of  the  i  ”  ply 
material  axes 

vector  of  thermal  coefficients  of  expansion  of  the  i*^*^  ply  referred  to  ply's 
material  axes 

ply  longitudinal  thermal  coefficient  of  expansion 
ply  transverse  thermal  coefficient  of  expansion 
correlation  coefficient  for  ply  transverse  tensile  strength 
void  strain  magnification  coefficient 

vector  of  composite  strains  referred  to  composite  structural  axes  at  the  ref¬ 
erence  plane 

vector  of  strains  for  the  i*"^  ply  referred  to  ply's  material  axes 
matrix  allowable  strain  in  tension 

ply  angle  measured  from  the  composite  structural  axes  to  the  ply  material 
axis  1 

composite  local  curvatures  referred  to  composite  structural  axes 

ply  Poisson's  ratio  in  the  two  direction  when  the  load  is  applied  in  the  one 
direction 

vector  of  ply  stresses  referred  to  ply's  material  axes 
strain  magnification  factor  for  ply  transverse  strength 
vector  or  column  matrix 
square  array 


GOVERNING  EQUATIONS 

The  laminate  theory  equation  for  predicting  the  ply  residual  stresses  is  given  by 

I'^il  1  =  [Eiif'dRiiHwl  -  <<cxi  -  (C;,  }> 

The  reference  plane  strains  and  the  curvature  changes  {k_^}  for  a  free  com- 

COX  cx 


5 


posite  (no  external  loads  and  no  boundary  constraints)  are  computed  from 


’[^cxlICcxl’ 

-1 

{ 

U^cx  '  J 

PcJlI>cxl. 

l^^cATx  U 

(2) 


Ecjuations  (1)  and  (2)  show  that  the  ply  residual  stress  is  a  function  of  the  following 
factors:  (1)  the  composite  structural  stiffness,  (2)  the  ply  spatial  location  in  the  com¬ 
posite,  (3)  the  ply  stress-strain  relations,  (4)  the  ply  thermal  coefficients  of  expansion, 
and  (5)  the  temperature  difference  between  the  ply  and  the  reference  value.  This  differ¬ 
ence  equals  the  ply  temperature  minus  the  cure  or  processing  temperature  in  computing 
residual  stresses.  The  factors  (1),  (3),  and  (4)  are  related  to  the  ba.sic  variables  (con¬ 
stituent  materials  and  thermoelastic  properties,  fiber  and  void  volume  ratios,  number 
of  plies,  and  ply  orientation).  These  relations  are  given  in  either  reference  10  or  14. 


RESIDUAL  STRESSES  IN  BORON/A LUMINUM  COMPOSITES 

The  temperature  difference  selected  for  boron/aluminum  composites  is  900^  F 
(773  K).  The  room  temperature  properties  are  used  in  the  linear  analysis.  The  appli¬ 
cability  of  linear  analysis  and  room-temperature  properties  might  be  questioned  for  this 
temperature  difference,  particularly  for  the  aluminum  rnalrix.  However,  there  are  four 
reasons  that  justify  this  approach.  First,  the  residual  stresses  are  computed  at  room 
temperature.  Second,  the  aluminum  matrix  modulus  decreases  with  increasing  temper¬ 
ature,  but  its  thermal  coefficient  of  expansion  increases.  These  two  effects  tend  to  be 
compensatory.  The  third  rea.son  is  thal  the  results  from  a.  linear  analysis  can  be  scaled 
down  proportionally  to  any  desired  temperature  difference.  The  fourth  reason  is  that 
composites  that  are  designed  to  carry  the  residual  stress  at  room  temperature  caused 
by  a  900^  F  (773  K)  temperature  difference  will  be  conservatively  designed. 

Ply  residual  transverse  and  shear  stresses  (as  defined  in  fig.  2)  are  plotted  as 
functions  of  ply  orientation  angle  and  fiber  volume  ratio  in  figure  3 .  The  composite 
geometry  and  the  temperature  difference  is  noted  in  the  captions  of  these  figures. 

Ply  residual  transverse  stresses  as  a  function  of  ply  angle  are  plotted  in  figure  4 
for  the  general  composite  geometry  8[0,  2(±0),  2(90),  2(+0),  0].  The  temperature  dif¬ 
ference  and  the  fiber  volume  ratio  (FWO  are  900^  F  (773  K)  and  0.  55,  respectively. 

The  results  in  this  figure  for  6  -  0^  or  90^  illustrate  the  dependence  of  the  ply  residual 
stress  on  the  number  of  plies  in  each  orientation. 

Ply  residual  transverse  stresses  for  composites  of  the  following  three  geometries 
8[4(+0),  4(-0)],  8(±0),  8[4(±0),  4(t0)]  and  with  an  FVR  of  0,  55  are  plotted  against  6  in 


6 


figure  5.  The  ply  considered  in  these  composites  is  the  first  ply  in  the  sequence.  The 
composite  geometry  8[4(+0),  4(-0)]  represents  a  composite  that  is  nonsymmetric  with 
respect  to  bending.  The  composite  geometry  8(±0)  represents  a  composite  that  is  non¬ 
symmetric  with  respect  to  bending  with  interspersed  plies.  The  composite  geometry 
8[4(±0),  4(+0)]  represents  a  composite  that  is  symmetric  with  respect  to  bending  with 
interspersed  plies. 

The  results  in  figures  3  to  5  show  that  the  ply  residual  transverse  and  shear 
stresses  increase  with  ply  angle  but  decrease  with  fiber  volume  ratio.  The  smaller  the 
number  of  plies  in  a  given  orientation,  relative  to  other  orientations  in  the  same  com¬ 
posite,  the  greater  the  residual  stresses  in  these  plies.  The  residual  stresses  in  the 
outer  plies  of  nonsymmetric  noninterspersed  composites  are  quite  high.  They  are  about 
twice  as  high  as  those  in  either  nonsymmetric  interspersed  or  symmetric  interspersed 
plies.  Introduction  of  ±45°  plies  for  increased  torsional  stiffness  produces  high  residual 
stresses.  The  residual  stresses  in  quasi (pseudo)-isotropic  symmetric  composites  (such 
as  8[0,  2(±45),  2(90),  2 (+45),  0],  fig.  6)  are  equal  in  all  the  plies. 

The  transverse  ply  tensile  strength  for  a  boron/aluminum  unidirectional  composite 
(ply)  at  an  FVR  of  0.  55  is  14. 1  ksi  (97. 1  MN/m^)  (see  table  III).  This  value  is  exceeded 
in  many  of  the  composites  shown  in  figures  3  to  5. 

To  what  extent  the  residual  stresses  can  be  annealed,  if  any,  and  how  much  these 
stresses  affect  the  in-situ  ply  strength  is  not  known.  Also,  the  effect  of  creep  in  re¬ 
lieving  residual  stresses  is  not  known  at  this  time. 


RESIDUAL  STRESSES  IN  BORON,  S-GLASS,  AND 
THORNE  L- 50  GRAPHITE /EPOXY  COMPOSITES 

The  residual  stresses  in  these  composites  are  computed  from  a  300°  F  (440  K)  tem¬ 
perature  difference,  room  temperature  constituent  material  properties,  and  the  linear 
theory  available  in  reference  14.  This  theory  is  briefly  described  in  the  INTRODUC¬ 
TION.  Five  reasons  justify  this  approach:  (1)  The  residual  stresses  are  computed  at 
room  temperature.  (2)  Polymerization  shrinkage  (matrix  phase  change),  though  not 
accounted  for  directly  in  the  theory,  is  compensated  for  by  the  higher  temperature. 

(3)  The  epoxy  modulus  decreases  with  increasing  temperature,  but  the  epoxy  thermal 
coefficient  of  expansion  increases.  (4)  The  residual  stress  results  presented,  herein, 
can  be  scaled  down  proportionally  to  any  desired  temperature  difference.  (5)  Com¬ 
posites  designed  to  withstand  the  residual  stress  at  room  temperature  predicted  by  this 
approach  will  be  conservatively  designed. 

Results  for  ply  residual  transverse  stresses  are  presented.  These  stresses  are 
primarily  responsible  for  the  transply  cracking. 
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Boron/Epoxy  Composites 


The  residual  transverse  stresses  in  the  ±6  plies  of  the  boron/epoxy  composites 
8[2(0),  2(±0),  2(+6),  2(0)]  are  plotted  against  8  and  three  fiber  volume  ratios  in  fig¬ 
ure  6(a). 

Ply  residual  transverse  stresses  for  the  0°,  ±8°,  and  90°  plies  of  the  composites 
8[0,  2(±e),  2(90),  2(+0),  0]  are  plotted  against  8  in  figure  6(b).  The  FVR  is  0.  55  for  these 
composites.  The  results  in  this  figure  at  0  =  0°  and  90°  illustrate  the  effect  on  re¬ 
sidual  stress  of  the  numbers  of  plies  in  each  orientation. 

The  residual  transverse  stress  in  the  outer  ply  of  the  composites  8[4(+0),  4(-0)], 
8(±0),  and  8[4(±0),  4(+0)]  is  plotted  against  0  in  figure  6(c).  The  FVR  is  0.  55  for  these 
composites.  The  composites  8[4(+0),  4(-0)]  are  nonsymmetric  with  respect  to  bending 
and  have  a  noninterspersed  ply  stacking  sequence.  The  composites  8(±0)  are  nonsym¬ 
metric  with  respect  to  bending  but  do  have  an  interspersed  ply  stacking  sequence.  The 
composites  8[4(±0),  4(+0)]  are  symmetric  and  have  interspersed  ply  sequence. 

The  transverse  strength  for  the  unidirectional  composite  (ply)  at  an  FVR  of  0.  55  is 
8.  6  ksi  (59.  3  MN/m^)  (see  table  III).  This  value  is  exceeded  only  by  the  residual 
stresses  in  the  outer  plies  of  the  composite  8[4(+0),  4(-0)]  (fig.  6(c)).  However,  the 
margin  of  safety  is  about  0.  3  in  many  of  the  other  composites. 

Several  interesting  points  are  illustrated  by  the  results  shown  in  figure  6:  (1)  The 
ply  residual  transverse  stress  in  boron/epoxy  composites  is  relatively  high  when  com¬ 
pared  with  the  ply  transverse  strength  for  some  composite  geometries.  (2)  The  ply  re¬ 
sidual  stress  is  quite  sensitive  to  the  ply  angle,  but  not  as  sensitive  to  FVR.  (3)  The 
ply  transverse  residual  stresses  in  the  composites  8[0,  2(±0),  2(90),  2(+0),  0]  are  practi¬ 
cally  invariant  with  0.  (4)  Nonsymmetric  noninterspersed-ply  composites  experience 

high  ply  residual  stresses  when  compared  with  other  composite  geometries.  (5)  Non¬ 
symmetric  interspersed-ply  and  symmetric  interspersed-ply  composites  experience  ply 
residual  stresses  of  about  the  same  magnitude. 


S- Glass/Epoxy  Composites 

The  residual  transverse  stress  in  the  ±0  ply  of  the  S-glass/epoxy  composites 
8[2(0),  2(±0),  2(+0),  2(0)]  is  plotted  against  0  and  three  FVR  in  figure  7(a). 

Ply  residual  transverse  stresses  for  the  0°,  ±0°,  and  90°  plies  of  the  composites 
8[0,  2(±0),  2(90),  2(+0),  0]  are  plotted  against  0  in  figure  7(b).  The  FVR  is  0.  55  for 
these  composites.  The  results  in  this  figure  for  0  =  0°  and  90°  illustrate  the  effect  on 
residual  stress  of  the  number  of  plies  in  each  orientation. 

The  residual  transverse  stresses  in  the  outer  ply  of  the  composites  8[4(+0),  4(-0)], 
8(±0),  and  8[4(±0),  4(+0)]  are  plotted  against  0  in  figure  7(c).  The  FVR  is  0.  55  for 
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these  composites.  The  composites  8[4(+e),  4(-0)]  are  nonsymmetric  noninterspersed- 
ply  composites.  The  composites  8(±0)  are  nonsymmetric  interspersed- ply  composites. 
The  composites  8[4(±e),  4(+0)]  are  symmetric  interspersed- ply  composites. 

The  transverse  strength  for  the  S-glass/epoxy  unidirectional  composite  (ply)  at  an 
FVR  of  0.  55  is  4.  4  ksi  (30.  4  MN/m^)  (see  table  III).  This  value  is  very  close  to  the  ply 
residual  stress  of  several  composites.  It  is  exceeded  in  the  outer  plies  of  some  com¬ 
posites  8[4(+0),  4(-0)].  This  implies  that  S-glass/epoxy  composites  of  the  geometries 
and  FVR  examined  in  figure  7  could  experience  transply  cracks.  This  will  not  be  the 
case  if  the  transverse  strength  of  the  plies  is  much  greater  than  4.  4  ksi  (see,  e.  g. , 
table  1  of  ref.  5). 

Several  interesting  points  are  illustrated  by  the  results  in  figure  7.  (1)  The  calcu¬ 

lated  ply  residual  transverse  stress  could  exceed  the  ply  transverse  strength.  The  ex¬ 
tent  that  the  ply  transverse  strength  is  exceeded  by  the  calculated  ply  residual  stress 
depends  on  composite  geometry.  (2)  The  ply  residual  stress  is  quite  sensitive  to  ply 
orientation  angle  in  certain  composite  geometries.  (3)  The  ply  residual  stress  is  rela¬ 
tively  insensitive  to  fiber  volume  ratio  in  composites  with  less  than  0.  55  FVR.  (4)  The 
ply  residual  stress  is  sensitive  to  the  number  of  plies  in  each  orientation.  (5)  Nonsym¬ 
metric  noninterspersed-ply  composites  experience  comparatively  large  ply  residual 
stresses.  (6)  Nonsymmetric  interspersed-ply  and  symmetric  interspersed-ply  com¬ 
posites  experience  ply  residual  stresses  of  about  the  same  magnitude. 


Thornel-50  Graphite/Epoxy  Composites 

The  residual  transverse  stress  in  the  ±6  ply  of  the  Thornel-50  graphite/epoxy 
composites  8[2(0),  2(±0),  2(+0),  2(0)]  is  plotted  against  B  and  three  values  of  FVR  in  fig¬ 
ure  8(a). 

Ply  residual  transverse  stresses  for  the  0®,  ±6°,  and  90°-plies  of  the  composites 
8[0,  2(±6>),  2(90),  2(-0),O]  are  plotted  against  9  in  figure  8(b).  The  FVR  is  0.  55  for  these 
composites.  The  results  in  this  figure  for  6=0°  and  90°  illustrate  the  effect  on  re¬ 
sidual  stress  of  the  number  of  plies  in  each  orientation. 

The  residual  transverse  stress  in  the  outer  ply  of  the  composites  8[4(+6),  4[-6)], 
8(±6),  and  8[4(±6),  4(+6)]  is  plotted  against  6  in  figure  8(c).  The  FVR  is  0.  55  for  these 
composites.  The  composites  8[4(+0),  4(-0)]  are  nonsymmetric  noninterspersed-ply 
composites.  The  composites  8(±0)  are  nonsymmetric  interspersed-ply  composites. 

The  composites  8[4(±0),  4(+6)]  are  symmetric  interspersed-ply  composites. 

The  transverse  strength  for  the  Thornel-50  graphite /epoxy  unidirectional  composite 
(ply)  at  an  FVR  of  0.  55  is  3.  7  ksi  (25.  5  MN/m^)  (see  table  III).  This  value  is  exceeded, 
or  nearly  so,  by  several  composites.  This  implies  that  plies  in  Thornel-50  graphite/ 


9 


epoxy  composites  of  the  geometries  and  values  of  FVR  examined  in  figures  8(a)  to  8(c) 
will  experience  transply  cracks.  This  has  been  observed  experimentally  (refs.  5  to  7). 
This  will  not  be  the  case  if  the  transverse  strength  of  the  plies  is  about  8.  0  ksi 
(55.  2  MN/m^). 

Several  interesting  points  are  illustrated  by  the  results  in  figure  8.  (1)  The  calcu¬ 

lated  ply  residual  transverse  stress  exceeds  the  ply  transverse  strength.  The  extent 
that  the  ply  transverse  strength  is  exceeded  by  the  calculated  ply  residual  stress  de¬ 
pends  on  composite  geometry.  (2)  The  ply  residual  stress  is  quite  sensitive  to  ply  angle 
in  certain  composite  geometries.  (3)  The  ply  residual  stress  is  sensitive  to  FVR. 

(4)  The  ply  residual  stress  is  insensitive  to  the  number  of  plies  in  each  orientation.  It 
is  relatively  insensitive  to  9  in  the  composites  8[0,  2(±e),  2(90),  2(+i9),  0].  (5)  Nonsym- 

metric  noninterspersed-ply  composites  experience  ply  residual  stress  that  is  about 
30  percent  higher  than  other  composite  geometries.  (6)  Nonsymmetric  interspersed-ply 
and  symmetric  interspersed-ply  composites  experience  ply  residual  stresses  of  about 
the  same  magnitude. 


SIMULATED  COMPOSITE  COMPONENTS 

The  residual  stresses  in  the  plies  of  two  simulated  composite  components  were 
studied.  The  composite  geometries  are  24[4(±30),  2(±15),  12(0),  2(+15),  4(+30)]  and 
24[4(±45),  2(±22.  5),  12(0),  2(T22.  5),  4(+45)].  The  composite  system  is  Modmor-I  graphite/ 
polyimide,  and  it  is  used  for  elevated  temperature  application.  The  temperature  differ¬ 
ence  for  this  system  is  about  600°  F  (606  K). 

The  dependence  of  the  ply  residual  stresses  on  several  factors  is  examined.  These 
factors  are  (1)  fiber  volume  ratio  (FVR),  (2)  void  content,  (3)  variation  of  the  FVR  from 
ply  to  ply  within  the  same  composite,  (4)  different  groups  of  plies  with  different  FVR, 

(5)  variations  in  the  matrix  thermal  coefficient  of  expansion,  and  (6)  variations  in  the 
matrix  modulus.  Also,  the  effect  of  the  residual  stresses  on  the  relative  rotation  of 
adjacent  plies  is  investigated. 

The  computations  are  based  on  a  600°  F  (606  K)  temperature  difference  and  room- 
temperature  constituent  material  properties.  The  justification  for  this  approach  is  the 
same  as  that  given  for  the  other  fiber /epoxy  composite  systems.  (See  section  Residual 
Stresses  in  Boron,  S- Glass,  and  Thornel- 50/Epoxy  Composites. ) 

The  ply  residual  transverse  stress  is  plotted  against  FVR  in  figure  9(a)  for  both 
components.  Results  for  the  ply  residual  longitudinal  stress  are  shown  in  figure  9(b). 

The  residual  longitudinal  stress  in  the  0°  plies  of  both  components  is  tensile.  Results 
for  ply  residual  shear  stress  are  plotted  against  FVR  in  figure  9(c).  The  results  in 
these  figures  show  that  the  ply  residual  stresses  are  sensitive  to  both  composite  geom¬ 
etry  and  FVR. 
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The  void  content  effect  on  the  ply  transverse  residual  stress  is  illustrated  in  fig¬ 
ure  10.  This  effect  is  negligible  for  all  practical  purposes.  However,  it  is  important 
to  keep  in  mind  that  the  ply  transverse  strength  decreases  with  increasing  void  content 
(ref.  12).  The  decrease  in  ply  transverse  strength  could  produce  transply  cracks  at 
smaller  temperature  differences.  Proper  consideration  of  the  void  effects  is  necessary 
to  avoid  transply  cracking.  (See  eq.  (3)  in  this  section. ) 

The  effects  of  variable  ply  FVR  on  the  ply  transverse  residual  stresses  are  shown 
in  figure  11  for  one  component  only.  These  effects  are  negligible  as  can  be  seen  in  the 
figure.  Having  different  groups  of  plies  at  different  FVR  has  a  negligible  effect  on  the 
ply  residual  transverse  stress.  This  is  illustrated  in  figure  12.  Negligible  effects  were 
foimd  also  when  the  ply  temperature  was  varied  similar  to  the  FVR.  These  results  are 
not  presented  here. 

The  dependence  of  the  ply  residual  transverse  stress  on  the  matrix  thermal  coeffi¬ 
cient  of  expansion  is  illustrated  in  figure  13.  Two  important  points  are  observed  here: 
(1)  The  ply  residual  stress  varies  linearly  with  matrix  thermal  coefficient  of  expansion 
and  (2)  the  residual  stresses  can  be  reduced  by  decreasing  the  magnitude  of  the  thermal 
coefficient  of  the  matrix.  The  effects  of  varying  the  fiber  thermal  coefficients  of  expan¬ 
sion  will  produce  an  effect  opposite  to  those  shown  in  figure  13  for  the  thermal  coeffi¬ 
cient  of  the  matrix.  This  implies  that  increasing  the  fiber  thermal  coefficient  of  expan¬ 
sion  is  another  effective  way  to  reduce  the  residual  stress. 

The  dependence  of  the  ply  residual  transverse  stress  on  the  matrix  modulus  is  illus¬ 
trated  in  figure  14.  The  results  in  this  figure  show  that  the  residual  stress  increases 
nonlinear ly  with  increasing  matrix  modulus.  This  increase  is  quite  rapid  at  modulus 
values  of  100  to  700  ksi  (690  MN/m^)  and  not  as  rapid  at  modulus  values  greater  than 
700  ksi  (4.  8  MN/m^). 

The  results  in  figure  14  should  be  examined  jointly  with  the  ply  transverse  strength 
before  any  meaningful  conclusions  are  made.  The  ply  transverse  strength  is  given  by 
the  expression  (ref.  12). 


-  h2'T~r~^^l22 
Pv‘^m22 

where  ^i22T  transverse  tensile  strength,  correlation  factor, 

e  is  the  allowable  matrix  tensile  strain  to  denote  failxire,  /3  is  the  void  strain 
magnification  factor,  (P ^22  strain  magnification  factor,  and  E^22 

transverse  modulus.  The  variables  ^22’  ®Z22  independent  of  the 

matrix  modulus.  For  example,  <P^j22  decreases  but  £^£2  increases  with  increasing 
matrix  modulus  (ref  10),  and  e^^prp  increases  with  increasing  matrix  modulus  in  gen¬ 
eral. 
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The  follov/ing  point  can  be  made  in  view  of  the  preceding  discussion.  The  ply  trans¬ 
verse  residual  stresses  increase  v/ith  increasing  matrix  modulus.  The  ply  transverse 
strength  also  increases  v/ith  increasing  matrix  m.odulus.  Transply  cracking  will  even¬ 
tually  occur  if  the  ply  residual  transverse  stress  increases  at  a  higher  rate  than  the  ply 
transverse  strength.  Limited  results  showed  thad  the  ply  transverse  strength  increases 
a,t  a  higher  rate  than  the  residual  stress  for  some  graphite  composites  during  this  inves¬ 
tigation. 

The  effect  of  the  residual  stress  on  the  adjacent  ply  relative  rotation  is  illustrated 
in  figure  15.  Relative  rotation  is  plotted  against  fiber  volume  ratio  for  selected  adjacent 
plies  of  the  tw'o  components.  The  estimated  allov/able  relative  rotation  is  also  plotted  in 
the  figure,  and  the  prediction  of  it  is  discussed  in  references  10  and  11.  The  results  in 
figure  15  shov/  that  residual  stresses  could  cause  interply  de laminations.  The  de lami¬ 
nation  depends  on  the  constituent  materials  amd  composite  geometry,  the  ply  angles  of 
the  two  adjacent  plies,  and  the  fiber  and  void  volume  ratios.  It  is  interesting  to  note 
that  the  component  with  the  ±30^  and  ±15^  plies  has  a  higher  probability  of  interply  de¬ 
lamination  than  the  component  with  the  ±45^  and  ±22,  5^  plies.  The  component  with  the 
±30*^  and  ±15^  plies  is  superior  in  ail  other  respects  with  respect  to  residual  stress.  To 
the  author^s  Imovdedge,  the  significance  of  the  results  in  figure  25  has  not  been  recog¬ 
nized  by  other  investigators  to  date. 

The  results  in  figures  16  to  25  show  that  com.posites  designed  by  neglecting  the  re¬ 
sidual  stresses  vdll  be  unconservaTive  in  some  loading  conditions. 

The  effects  of  introducing  transitional  plies  (intermediate  plies  for  less -abrupt 
orientation  transition  betv/een  groups  of  plies)  or  decreasing  the  number  of  off-axis  plies 
on  the  residual  stress  are  illustrated  in  figure  18.  The  residual  transverse  stress  in 
the  ±45^  plies  is  plotted  agadnst  FVR.  The  composite  with  the  transitional  plies  is  de¬ 
noted  by  the  long-dashed  curve.  The  composite  in  v/hich  the  transitional  plies  have  been 
incorporated  v/ith  the  face  (surface)  plies  is  denoted  by  the  solid  curve.  The  composite 
in  v/hich  the  transitionad  plies  have  been  incorporated  in  the  core  plies  is  identified  by 
the  short-dashed  curve. 

The  results  in  figure  16  show  that  both  reduction  of  the  number  of  off-axis  plies  and 
introduction  of  transitional  plies  reduce  the  residual  stress  by  small  amounts.  The 
number  of  off-axis  plies  is  controlled  by  strength  and  stiffness  requirements  and  changes 
might  not  be  permitted  in  specific  designSc 


WAYS  TO  MINIMIZE  OR  ELIMINATE  LAMINATION  RESIDUAL  STRESSES 


Several  alternatives  are  possible  for  minimizing  or  elimdnating  the  lamination  re¬ 
sidual  stresses.  These  alternadives  may  be  grouped  into  two  categories.  The  first 
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category  is  general,  and  the  second  is  for  specific  designs.  They  are  summarized 
briefly. 

In  general,  residual  stresses  at  room  temperature  can  be  reduced  or  eliminated  by 

(1)  Use  of  lower  cure  or  processing  temperature  matrices 

(2)  Decreasing  the  difference  between  the  fiber  and  the  matrix  thermal  coefficients 

of  expansion 

(3)  Increasing  the  allowable  elongation  of  the  matrix  at  failure  with  a  simultaneous 

increase  in  the  matrix  modulus . 

In  specific  designs,  the  residual  stress  can  be  reduced  by 

(1)  Proper  selection  of  composite  system 

(2)  Suitable  ply  orientation  and  stacking  sequence 

(3)  Appropriate  selection  of  the  number  of  plies  in  each  direction 

(4)  Introduction  of  transitional  plies 

(5)  Decreasing  the  number  of  off-axis  pUes  if  stiffness  and  strength  requirements 

permit  it 

(6)  Increasing  the  fiber  volume  ratio 

(7)  Any  possible  combination  of  these. 


EXPERIMENTAL  WORK  NEEDED 

Both  theory  and  experiment  show  that  the  ply  residual  transverse  stresses  due  to 
lamination  are  of  magnitudes  comparable  to  those  of  the  ply  transverse  strengths.  How¬ 
ever,  four  important  questions  need  to  be  answered  before  the  effect  of  residual 
stresses  on  the  composite  structural  response  can  be  realistically  assessed:  (1)  decay 
of  residual  stress  with  time,  (2)  effect  of  transply  cracks  on  composite  stiffness, 

(3)  effect  of  transply  cracks  on  composite  strengths,  and  (4)  effect  of  transply  cracks  on 
composite  degradation  by  the  environment. 

These  questions  can  be  answered  best  by  simple  but  strategic  experimental  investi¬ 
gation.  Bending  tests  in  fatigue  and  creep  should  provide  qualitative  answers  to  these 
questions. 


CONCLUSIONS 

The  results  of  this  theoretical  investigation,  which  is  based  on  linear  laminate 
theory  and  on  a  semiempirical  micromechanics  theory,  lead  to  the  following  conclu¬ 
sions. 

1.  The  ply  residual  transverse  stress  is  generally  tensile.  The  residual  stress 
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depends  on  the  composite  system,  ply  orientation,  fiber  volume  ratio,  and  processing 
conditions.  The  residual  stress  can  be  large  relative  to  transverse  strength  for 
graphite /epoxy  and  boron/aluminum  composites. 

2.  The  ply  residual  in-plane  shear  stress  can  be  comparable  to  the  ply’s  shear 
strength. 

3.  The  ply  residual  transverse  or  shear  stress  can  be  comparable  to  the  corre¬ 
sponding  ply  strengths  in  simulated  components.  This  can  cause  transply  cracking.  It 
can  also  reduce  the  ply  capacity  to  carry  transverse  tensile  or  in-plane  shear  loads. 

4.  The  ply  residual  stresses  are  very  sensitive  to  temperature  difference,  ply 
orientation,  number  of  plies  in  each  direction  and  to  fiber  volume  ratio. 

5.  The  ply  residual  stresses  are  relatively  insensitive  to  void  content  and  to  mod¬ 
erate  ply-to-ply  fiber  volume  ratio  and  processing  temperature  variations. 

6.  The  ply  residual  stresses  generally  increase  with  increasing  ply  orientation 
angle . 

7.  The  ply  residual  stresses  increase  nonlinearly  with  increasing  matrix  modulus. 
They  increase  linearly  with  the  matrix  coefficients  of  expansion. 

8.  It  is  possible  to  orient  plies  in  simulated  components  so  that  some  plies  have 
tensile  longitudinal  residual  stress.  However,  the  residual  longitudinal  stress  is 
usually  compressive.  Its  magnitude  is  usually  15  to  20  percent  of  the  ply’s  correspond¬ 
ing  compressive  strength. 

9.  Transply  residual  stress  cracking,  in  general,  can  be  prevented  by  increasing 
the  ply  transverse  and  shear  strengths,  and  decreasing  the  processing  temperature.  It 
can  be  prevented  by  decreasing  the  difference  between  the  thermal  coefficients  of  the 
fiber  and  the  matrix  sufficiently. 

10.  In  specific  designs,  cracking  can  be  prevented  by  suitable  ply  orientation,  var¬ 
iation  of  the  number  of  plies  in  each  direction,  introduction  of  transitional  plies,  reduc¬ 
tion  of  the  number  of  off-axis  plies,  increasing  the  fiber  content  (within  practical  limits), 
or  by  any  combination  of  these. 

11.  Neglecting  the  residual  stresses  in  designing  with  fiber  composites  can  lead  to 
unconservative  design. 

12.  Experimental  work  is  needed  for  a  realistic  assessment  of  residual  stresses  on 
composite  structural  response. 

Lewis  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Cleveland,  Ohio,  October  5,  1970, 

129-03. 
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Figure  1.  -  Typical  fiber  composite  geometry  (plane  ABCD  is  the  reference  plane) 


Figure  2.  -  Single  ply. 


Fiber  volume 
ratio 


residual  stress  for  boron/aluminum  composite.  Composite  geometry,  8[2(0),  2(i0),  2(+0),  2(0)];  temperature  difference, 


Piy  orientation  angle,  6,  deg 

Figure  4.  -  Piy  residua!  transverse  stress  for  boron  aluminum  composites.  Com¬ 
posite  geometry,  8[0,  2(:i:8).  2(90),  2{~“8),  O] ;  fiber  volume  ratio,  0. 55;  tempera¬ 
ture  difference,  930°  F(773  K). 


Ply  orientation  angle,  6,  deg 

Figure  5.  -  Outer-ply  residual  transverse  stress  for  boron  aluminum  composites. 
Fiber  volume  ratio,  0.55;  temperature  difference,  900'  F  (773  K), 
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Composite  geometry 


Ply  orientation  angle,  0,  deg 
(c)  Outer  ply;  fiber  volume  ratio,  0.55. 

Ply  transverse  residual  stress  for  boron/ epoxy  composites.  Temperature  difference,  300“  F  (440  K). 


Fiber  volume 
ratio 


(c)  Outer  ply;  fiber  volume  ratio,  0.55. 

residual  transverse  stress  for  graphite  Thornel •■50/epoxy  composites.  Temperature  difference, 


Stress,  MN/m^ 


orientation 

angle, 

0. 

deg 


Fiber  volume  ratio,  kj 


(a)  Transverse  stress. 


Fiber  volume  ratio,  kf 


(b)  Shear  stress.  (q)  Longitudinal  stress. 

Figure  9.  -  Ply  residual  stress  for  Modmor-I/polyimide  composites.  Temperature  difference,  600°  F  (606  K) 


Figure  12.  -  Effect  of  having  different  groups  of  plies  at  different  fiber 
contents  on  ply  transverse  residual  stress.  Modmor-I/polyimide 
composite;  composite  geometry,  24  [4(±45),  2(±22.5),  12(0),  2{+22.5), 
4(1^45)];  zero  void  content;  temperature  difference,  600°  F  (606  K). 


Composite  geometry 
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- 24[4(±45),  2(±22.5),  12(0),  2(^22.5)41^45)] 

100  15  I — 


Ply 

orientation 


.  ^15 

v_30 


0l - ^ ^ - 1 

10  20  30  40x10  ° 

Matrix  thermal  coefficient  of  expansion,  in./in./°F 


Matrix  thermal  coefficient  of  expansion,  cm/cm/K 


Figure  13.  -  Effects  of  matrix  thermal  coefficient  of  expansion  on  ply  transverse 
residual  stress.  Modmor-l/polyimide  composites;  zero  void  content;  fiber 
volume  ratio,  0.50;  temperature  difference,  600°  F  (606  K). 


Relative  rotation,  deg 


Composite  geometry 

- 24[6(±45),  12(0},  6(t45)] 

" - 24[4{±45),  2(±22.5),  12(0},  2(+22.5),  4(+45)] 

- 24[4(±45),  16(0),  4(T45)] 


.4  .6  .8 

Fiber  volume  ratio, 

gure  16.  "  Transitional  plies  and  decreasing  number  of  off-axis 
plies  effect  on  J:45''  ply  residual  transverse  stress.  Modmor-I/ 
polyimide  composite;  temperature  difference,  600"'  F  (606  K). 
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